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Detection of electronic nematicity using scanning tunneling microscopy 
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Electronic nematic phases have been proposed to occur in various correlated electron systems and were 
recently claimed to have been detected in scanning tunneling microscopy (STM) conductance maps of 
the pseudogap states of the cuprate high-temperature superconductor Bi2Sr2CaCu208+(5 (Bi-2212). We 
investigate the influence of anisotropic STM tip structures on such measurements and establish, with a 
model calculation, the presence of a tunneling interference effect within an STM junction that induces 
energy-dependent symmetry -breaking features in the conductance maps. We experimentally confirm this 
phenomenon on different correlated electron systems, including measurements in the pseudogap state of 
Bi-2212, showing that the apparent nematic behavior of the imaged crystal lattice is likely not due to 
nematic order but is related to how a realistic STM tip probes the band structure of a material. We further 
establish that this interference effect can be used as a sensitive probe of changes in the momentum structure 
of the sample's quasiparticles as a function of energy. 

PACS numbers: 74.72.Kf, 74.70.Tx, 71.27.+a, 74.55.+V 
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The concept of broken symmetry is essential to con- 
densed matter physics. Identification of the fundamental 
symmetries of a solid-state system leads to the understand- 
ing of the low-energy excitations which govern its prop- 
erties. For example, the unraveling of the three-decade old 
mystery of unconventional superconductivity hinges on de- 
termining the symmetries of the correlated electronic state 
from which cooper pairs are formed. Recently, electronic 
nematic phases, where electronic states undergo a spon- 
taneous four-fold (C4) to two-fold (C2) symmetry break- 
ing, have gained much interest as a possible candidate for 
various hidden order states in several correlated electron 
systems such as cuprates, iron-based superconductors, and 
heavy fermion materials |[I]-[6l. However, such states are 
difficult to detect using non-local probes because of pos- 
sible twin-domain structures in macroscopic samples. Re- 
cently STM has been proposed as the method of choice for 
the detection of four- fold electronic symmetry breaking fT"- 
[14|. Lawler et al [8|, and Messaros et al L9J in a subse- 
quent study, interpreted the rotational symmetry breaking 
in the STM data as evidence for an electronic nematic state 
inside the pseudogap phase of Bi-2212. 

The naive expectation has been that the influence of 
the STM tip's geometric structure is limited to induc- 
ing an easy-to-detect anisotropy in STM topographs or 
to influence energy-resolved STM differential conductance 
(dl/dV) maps in an energy independent manner. Here we 
show, through model calculations and experimental mea- 
surements on three correlated electron materials (CeCoIns, 
Bi-2212, and URu2Si2), that a tunneling interference ef- 
fect within an STM junction composed of a realistic tip 
(with some spatial anisotropy) can result in an artificial 
energy-dependent symmetry breaking of the STM conduc- 
tance maps. This phenomenon can occur even when the 



STM topograph taken with the same tip appears to be sym- 
metric. We demonstrate that previously reported two-fold 
symmetric conductance maps in high-Tc cuprates [E','9] are 
not evidence for rotational symmetry breaking (C4 to C2) 
originating from a nematic phase in these materials, but 
are rather due to the interference effect we have uncovered 
here. In this system, systematic measurements with dif- 
ferent tips on the same area of the sample, reported here, 
are also used to clearly demonstrate the lack of nematic 
order, without relying on any pseudogap-specific assump- 
tions about the tunneling process. We further show that 
the interference effect within the STM junction can never- 
theless be used as a sensitive tool to detect changes in the 
quasiparticle band structure as a function of energy. 

We start our discussion by considering how STM probes 
the electronic structure of a sample's surface. Following 
Tersoff and Hamann ifTSl the sample wave function can be 
written as 

^s,k(^) = ^a^exp[iKG ' r- (/^^ + I^gD^^] (1) 



where Rq — k ^ G defines the Wannier states, while 
Kj — \j2vn(\)lh is related to the work function for elec- 
tron decay into the vacuum, and the summation is over the 
reciprocal lattice vectors G [ 16|. Most discussions of STM 
data assume a metalic tip (energy-independent density of 
states), and approximate the STM differential conductance 
dl /dV (at small bias) as the spatial convolution (*) of the 
tip {pt) and sample {ps) densities of states (DOS) [ilSJ : 

-^{eV,r)(xpt^Ps{eV) (2) 

Wiihps{eV) = Y.S \^s\^5{Es-Ef-eV). Thenunderthe 
usual assumption of an isotropic tip, the conductance maps 
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FIG. 1. (Color online) (a) Parabolic band structure (E = 
^0 + tk'^, with t = (SO/tt^) and Eq = 0.59t) used to gener- 
ate (b), (c) and (d). (b) Simulated topography siieV = E = 3.95t 
imaged by the tip in the lower inset. The inset of the tip repre- 
sents the simulated |V^^p and is plotted on the same spatial scale 
as the lattice, (c) Simulated differential conductance at ^ = 3.8t 
showing rotational symmetry breaking with the same tip as in (b). 
(d) Simulated differential conductance at E = 2.3t showing no 
rotational symmetry breaking with the same tip as in (b). (e-g) 
represent the intensity (normalized to the maximum) of the two 
orthogonal Bragg peaks generated by the DFTs of (b), (c), and 
(d), respectively, (h) One-dimensional schematic representation 
of the interference between the wavefunction of a double tip with 
a quasiparticle state of momentum k. 



are simply proportional to the local DOS of the sample. 
Within this model of STM measurements the topograph is 
constructed from integrating such maps between the Fermi 
level up to the tip-sample bias. 

Before considering an anisotropic tip, we emphasize that 
Eq. ^ above is an approximation of Bardeen's formula for 
tunneling [[BUTTll 



Y^f{Es)[l - f{E, + eV)]\Ms,\'6{AE) (3) 



s,t 



where f{E) is the Fermi-Dirac distribution, and AE = 
Es — Ef. Here Mgt is the matrix element between the tip 



and the sample with the following spatial structure fTJ]: 

^^j d^-{rs^^,-^,vrs) (4) 



M, 



We now demonstrate the effect of the tip geometry on 
STM measurements, by simulating the conductance maps 
using Bardeen's matrix element (Eq. (|4]) above) with a two- 
fold symmetric tip structure, characterized by orthogonal 
lengths 5x = 0.2, 5?/ = 0.9 (lattice constant set to unity) 
|[T6l . We assume that the sample has a C4 symmetric 
electronic structure with a generic parabolic band structure 
Fig.fTla), which is isotropic in the k^ - ky plane. Calcu- 
lations of the STM topographic image (at eV — 3.95t) 
and conductance maps at two different energies of this 
four-fold symmetric sample with elongated tip wave func- 
tion are shown in Fig.[TJb-d). While the STM conduc- 
tance maps can show apparent asymmetry in the x - ^ plane 
Fig.[TJc), the topograph appears to be remarkably four-fold 
symmetric Fig.[T];b). Clearly, the summation of dl/dV 
maps over an appropriate range of energies can lead to a 
four-fold symmetric topograph. This finding demonstrates 
that a four- fold symmetric topograph cannot be used as an 
accurate method to characterize the STM tip geometry, as 
it is often assumed (for example see Ref. [8,]). 

An anisotropic tip would naturally induce an apparent 
breaking of four-fold symmetry in measurements of the 
electronic structure of a four- fold symmetric sample, as the 
conductance maps demonstrate. However, the energy de- 
pendence of the x-y asymmetry (which can change sign, 
see below), or its absence for some energies (Fig.[TJd)), 
points to a previously overlooked interference effect of 
STM measurements. Examining Eq. ^ we realize that the 
periodic corrugation along the x and y directions in the 
conductance maps are determined by the interference (see 
Fig.[TJh) for schematic) between the sample's quasiparticle 
states ^^ ^ with momentum k together with those of the 
tip (characterized in our two-fold symmetric tip by 6x and 
6y). Previous studies of the influence of asymmetric tips 
|[T8l - [20 l have overlooked the interference within the STM 
junction by using the approximation in Eo. ([2]) which ig- 
nores the phase information (e*^"^, see Eq. (jj)) of the sam- 
ple wave functions that are relevant in the evaluation of 
Eq. (|4]). Additionally, studies of quantum interference ef- 
fects in tunneling junctions [[2TI - [24ll have not considered 
the effects of geometrically asymmetric tips on the mea- 
surement of long-range periodic structures by the STM. 
In contrast, our model calculations clearly show that the 
electronic structure of a four- fold symmetric square lattice 
probed by a real STM tip can be two-fold symmetric de- 
pending on the momentum k (consequently energy) of the 
quasiparticles probed. 

For a more detailed analysis of the energy-dependence 
of this asymmetry in STM conductance maps, and to make 
a connection to experimental measurements, we quan- 
tify the calculated STM conductance maps with the two- 
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FIG. 2. (Color online) (a) DFT of dI/dV{E = 0.79t) generated 
using Eq. ([2| showing strong peaks due to the long-range periodic 
lattice strucutre. (b) Energy dependence of the lattice asymmetry 
parameter calculated from Eq. ^ for different tip configurations 
{5y = 0.3, and 5x as indicated in the figure). Notice that for a 
four-fold symmetric tip {5x — 0.3 curve in (b)) On{E) = for 
all energies. 
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dimensional asymmetry parameter that is commonly used 
in the context of nematic ordering |[3: 



0^{E) 



X{E) - Y{E) 
X{E) + Y{E) 



(5) 



where X{E) and Y{E) are the energy dependent ampli- 
tudes of the two Bragg peaks along the orthogonal direc- 
tions obtained from discrete Fourier transformation (DFT), 
as indicated in Fig.|2ja). A map with On — corresponds 
to a four-fold symmetric image, whereas Oat = 1 indi- 
cates an image with zero corrugation along either the x or 
y direction (as expected for example for one-dimensional 
stripes). Figure[2];b) shows On{E) for different tip config- 
urations. Despite the simplicity of the model band struc- 
ture, On{E) shows a significant energy dependence over 
the entire bandwidth and even a sign change. This illus- 
trates the sensitivity of this tip-induced interference effect 
to the band structure. Although the magnitude of On{E) 
can only be understood by a detailed knowledge of ^^^ its 
energy dependence acts as a detector of changes in the mo- 
mentum structure of the quasiparticle states. 

Our model calculation suggests that materials with 
changes in their electronic band structure as a function of 
energy (such as a rapid change of band dispersion) are 
likely to be good candidates for exhibiting the interfer- 
ence effect associated with asymmetric tips. A good ma- 
terial candidate for such a study is the heavy fermion com- 
pound CeCoIns, which crystalizes in the tetragonal crystal 
structure, ensuring the four- fold symmetry of its electronic 
states. Recent STM studies on CeCoIns have demonstrated 
that the electronic structure of this compound exhibits the 
development of a hybridization gap and associated heavy 
bands near the Fermi energy at low temperatures [25 1. Fig- 
ure|3] shows a topograph (a) of CeCoIns, a DFT of a con- 
ductance map on the same area (b) together with the STM 
spectrum as a function of energy (c), which demonstrates 
the presence of a hybridization gap in this compound near 



FIG. 3. (Color online) (a) Topograph of CeCoIns (setpoint con- 
dition at — 200mV and 1.6 nA) showing a square lattice. For en- 
hanced contrast the derivative of the data is shown, (b) DFT of 
conductance (dl/dV) map taken at —52 meV over the same FOV 
as (a) showing strong Bragg peaks representing the square lat- 
tice. Inset shows an enlargement of the bottom right Bragg peak, 
(c) Tunneling spectrum averaged over the area in (a), (d) Energy 
dependence of the Bragg peak intensity obtained via the DFT op- 
eration, (e) Asymmetry parameter calculated via Eq. ([5]). Inset 
of (e) represents the asymmtery parameter of the topograph ac- 
quired simultaneously to the conductance maps. 



the Fermi energy. Also shown in Fig.[3]is the intensity of 
the Bragg peaks in the conductance maps as a function of 
energy (d) and the asymmetry parameter On{E) (e) intro- 
duced above. 

Approaching the energy window near the Fermi level, 
where we expect strong changes in the band structure of 
CeCoIns due to hybridization of spd- and /-like electrons, 
we find a strong C4 -symmetry breaking of the STM con- 
ductance maps (Fig.|3fd,e)). Both the Bragg peaks and 
On{E) show an energy-dependent asymmetric behavior 
in tandem with the features of the STM spectrum. The ap- 
parent breaking of C4 symmetry in this experiment is not 
associated with nematic order in CeCoIns but rather probes 
the strong momentum dependence of the band structure of 
this compound near the Fermi level. In fact, it is remark- 
able that On{E) is sensitive to the most subtle features 
in the spectra as a function of energy (see dashed lines in 
Fig.[3jc)), which are associated with changes in the elec- 
tronic momentum structure as previously detected in the 
quasiparticle interference of this compound ll25l . 

Further evidence that the asymmetry between X and Y 
detected in the conductance maps of CeCoIns is associ- 
ated with interference in the STM junction can be found by 
repeating the same experiment with slightly different tips 
(created by interacting with the surface) over the same field 
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FIG. 4. (Color online) (a) On{E) measured on CeCoIns with different tips at 20 K. (b) On{E) measured on Bi2Sr2(Ca,Dy)Cu208+(5 
(Tc = 35 K) at 30 K over the same FOV with different tip configurations, (c) Open symbols represent On{E) measured on 
URu2Si2 above Tho (20 K) and below (15 K) over the same FOV with the same tip. Measurements on a second FOV with differ- 
ent tip configurations (closed symbols) were carried out below Tjjo (15 K). For comparison purposes, the average tunneling spectra 
are displayed (dashed curves) for the respective materials (for URu2Si2 the average spectrum at 10 K is displayed) |16|. The insets 
represent the asynmitery parameter of the respective individual topographs acquired simultaneously to the conductance maps. 



of view (FOV) or equivalent areas of the same cleaved sam- 
ple. As expected from our model calculations (Fig.[2];b)) 
different tips have different sensitivity to the momentum 
structure of the electronic structure of the sample, and de- 
pending on their geometry exhibit different degrees of C4- 
symmetry breaking in the acquired conductance maps. As 
Fig.Qa) shows, the energy dependence of the asymme- 
try parameter in the conductance maps, On{E), is a very 
strong function of the tip and not always correlated with 
the presence, or the degree of, Bragg peak asymmetry in 
the STM topographs of the same area. 

We turn our attention next to the claims that STM mea- 
surements of underdoped Bi-2212 samples break C4 sym- 
metry and exhibit nematic order mm. As shown in 
Fig.Qb) measurements on such a sample exhibit very sim- 
ilar characteristics to those of CeCoIns, where changes in 
the spectrum associated with the pseudogap coincide with 
apparent asymmetry and a non-zero On{E) in this energy 
window. Not only this correlation is very characteristic of 
the tip-induced symmetry breaking originating from inter- 
ference effects within the STM junction, we also find that 
On{E) displays a strong sensitivity to the tip structure 
when probing the exact same FOV with slightly different 
tips. Remarkably, tips that show very similar, nearly x-y 
symmetric, topographs can show very different energy de- 
pendences for On{E), and even exhibit opposite signs for 
the effect on the same exact area of the sample. Clearly, 
such behavior is more consistent with the tip-dependent in- 
terference in the STM junction, associated with changes 
in the momenta of electronic states within the pseudogap, 
rather than any nematic order. Consistent with this view. 



and with previous experiments lElIll, no domain bound- 
aries between regions of different nematic order parameter 
have ever been found despite the large areas used for STM 
studies. 

Before we conclude, we present experiments on one 
more materials system, the results of which demonstrate 
that the interference within the STM junction and the as- 
sociated asymmetry parameter can in fact be used to probe 
the onset of sudden changes in electronic band structures 
of materials. We have carried out temperature dependent 
experiments on the heavy fermion URu2Si2, which shows 
a sharp second order phase transition in the so-called "hid- 
den order state" below Tho — 17.5 K, the nature of which 
continues to be a mystery |[26l - [28ll . Experiments on this 
compound are also consistent with the asymmetry parame- 
ter picking up changes in the electronic states at low tem- 
peratures through the tip-dependent interference. However, 
contrasting measurements below and above the Tho, over 
the same FOV, and with the same tip (open symbols in 
Fig.Qc)), shows that the signals in On{E) change from 
a peak-like shape to a smooth curve, directly reflecting the 
change in the band structure as the hidden order phase tran- 
sition is crossed. At temperatures just below the transition, 
when the changes in the electronic states are difficult to de- 
tect in the STM spectra, we find that On{E) is extremely 
sensitive to the changes that occur in the electronic struc- 
ture of this material below Tho- 

Overall, our systematic measurements on three differ- 
ent materials (with three different characteristic gap en- 
ergy scales, 30 meV for CeCoIns, 4 meV for URu2Si2, and 
100 meV for Bi-2212), demonstrate the strong sensitivity 



of the asymmetry parameter On{E) to different tip con- 
figurations, and, specifically, how it can change sign for 
measurements over the same FOV. From these results, we 
can only conclude that On (E) is not a measure of the sym- 
metry breaking of the electronic states of the sample, rather 
it is the result of the interference effect which is evident in 
the elementary model of tunneling from a realistic tip dis- 
cussed in here. Although STM can in principle detect the 
onset of nematic order, we have demonstrated that sym- 
metry analyses of conductance maps can be dominated by 
the energy-dependence of the band structure of the sam- 
ple rather than nematic order. Perhaps the only experi- 
mentally reliable approach to detect rotational symmetry 
breaking order with the STM would be to image the pres- 
ence of domain boundaries, such as the structural ordering 
[10| and electronic smectic ordering [9| in Bi-2212, and 
electronic nematic ordering in iron-based superconductors 
171 HH. Alternatively rotation of the STM tip by an appro- 
priate angle (90° in the case of C^ symmetry) while main- 
taining the same location on the sample could be devel- 
oped to discount the role of the tip geometry. Regardless, 
the interference within the STM junction with realistically 
anisotropic tips described here shows that such measure- 
ments are a sensitive probe of the changes in the quasipar- 
ticle states of the sample as a function of energy, even when 
such changes might not be apparent in STM spectra. 
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THE BARDEEN MATRIX ELEMENT 

In Bardeen's formalism [1, 2] the tunnehng current is 
determined by 



27re 



E/(^^)[l - /(^^ + eV)]\Mst\'5{Es - E, 



s,t 



(1) 

where s(t) indexes the sample (tip), f{E) is the Fermi- 
Dirac distribution, and Mgt is a matrix element between 
the sample and tip states. The spatial structure of the STM 



image is therefore contained in M^^, which Bardeen [1] 
showed to be 



M. 
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where S is any closed surface in the vacuum separating the 
tip and the sample. We choose the surface of the sample to 
be the z = plane. Then the tip is above at coordinates 
(xt, yti^t)^ and S can be the < z < Zt plane. Manip- 
ulation of Eq. (2) shows that Mst resembles a probability 
current through the z surface: 



Mst{xt,yt,zt) 




^l{x,y,z^ 



diptix - xt,y -yt,z - Zt) 

dz 



-^l)t{x-xt,y-yt,z-zt) 



diljl{x,y,z) 
dz 



dx dy (3) 
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It can now be seen explicitily that the matrix element is 
dependent on the position of the tip. Equation (3) describes 
the spatial structure of the STM measurement and it can be 
written in a simplified way: 



Mst(xt,yt,zt) 



U: 



d^t 



^t 



d^s 



(4) 



2m \ ^ dz dz 

where the * operation denotes convolution in the x and y 
coordinates, and Mst{xt^ yt^Zf) is sl function of the tip co- 
ordinates. 



SAMPLE AND TIP WAVE FUNCTIONS 



Sample 

According to Tersoff and Hamann [2] the sample wave 
function is written as 

^s,ki^) = ^a^exp[iKG • r- {n^ + \Rg?)^ A (5) 



G 



i^G — k + G 
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(6) 

(7) 



The z-dependence depends on the work function (j) and on 
G. For large G the penetration of the corresponding Bloch 



waves into the vacuum decreases exponentially. Therefore 
we restrict the sum in Eq. (5) to the three highest order 
terms: G = 0, G = Gx, and G — Gy, where with- 
out loss of generality we have taken the simplest case of 
a square lattice defined by reciprocal lattice vectors G^ 
and Gy. These three highest order terms also represent the 
simplest wave function which can produce a square lattice 
in Eq. (5). Here we have taken the values a^ — 1, and 
^G ~ ^G ~ 2e*5. Though these coefficients determine 
the spatial phase of the periodic lattice structure, different 
relative values have been tested and do affect our conclu- 
sions. The results of the present study can also be easily 
expanded to other Bravais lattices, or other long-range pe- 
riodic structures, by appropriate choice of the reciprocal 
space vector basis. 



Tip 



First we consider an ideal tip, where the characteristic 
scale of its wave function is much smaller than the pe- 
riod of ^s, then tl^t can be modeled as an infinetely sharp 
delta function, and Eq. (5), together with Eq. (4), yield 
Mst ^ l^sP and consequently dl/dV oc Ps . This is 
the usual STM assumption that the diferential conductance 



is proportional to the local density of states of the sam- 
ple. However, Tersoff and Hamann [2] showed that an in- 
finetely sharp tip is not necessary to obtain atomic resolu- 
tion with the STM. We extend their results to accomodate 
for anisotropic tips. Considering the typical separation be- 
tween sample and tip (^ 5 A) we model the tip as an s- 
wave: 



i't 



^ia.^-T{x'+V^ + {z-ztY)2 



(8) 



where a is the phase of this wave function and V deter- 
mines its penetration into the vacuum [3]. We find that the 
simplest choice to model an anisotropic tip is to write it's 
wave function as a combination of four s-waves with their 
centers separated by 2 x 5x, and 2 x 6y: 



^t 
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We have chosen this particular functional form for ^pt 
because of its simplicity, though any other choice of 
anisotropic wave function in the x-y plane suffices to pro- 
duce the apparent rotational symmetry breaking described 
in the text. 



CALCULATION STRATEGY 

The caculation is performed in the following steps. First, 
for each constant energy contour of the band structure in 
Fig. 1, the momentum states k are represented by a circle 
in the k^-ky plane. We then slice that circle in A^ points 
equally spaced by angle, which are represented by unique 
values of k^ and ky. For each pair we calculate Mgt ac- 
cording to Eq. (4). Finally we add all the IM^^P (where 
the s states in the summation of Eq. (1) are the several i/j^ ^ 
states for a particular energy) to obtain the conductance 
map at the energy of interest. We repeat this procedure 
for each energy and construct the curves in Fig. 2(a) of the 
main text. 



ANGLE BETWEEN TIP ANISOTROPY AXIS AND THE 
LATTICE 

In Eq. 9 we implicitely assume that the anisotropy axis 
of the tip wavefunction is aligned with one of the crys- 
tal axes. This is rarely the case in an actual experiement. 
However, for most values of relative angles between the 
tip anisotropy axis and the crytal lattice, the induced x-y 
asymmetry in the conductance maps should still occur. We 
confirm this expectation by performing our calculation for 
different relative angles (see Fig. 1). Notice that for a value 
of 45° the tip geometry is unable to distinguish between 
the X and y axes of the sample, and therefore the value of 
On{E) is identically zero for all energies. 





1 2 3 
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FIG. 1. On{E) (left panel) for select angles between the crystal 
lattice (vertical arrow, right panel) and the tip anisotropy axis. 
For all curves the tip wave function on the right panel (shown in 
the length scale of a 2 lattice constants) was used with different 
relative angle to generate On (E) . 



TIP-CHANGING PROCEDURE ON CeCoIns AND 
URu2Si2 

The tip is offset from the measurement area and the con- 
figuration of its apex is modified by moving the tip toward 
the surface until we see a saturation in the tunneling cur- 
rent (beyond the limits of our pre-amplifier). This we call 
a controlled poke onto the surface. We have determined 
the optimal voltage applied to the scanner piezo empiri- 
cally. Because of the metallic nature of the CeCoIns and 
URu2Si2 surfaces, this controlled poke does not compro- 
mise the stability of the tunneling junction, nor does it 
modify the poking site beyond a few tens of angstroms. 
The change in the tip apex configuration is confirmed by 
subsequent topographic measurements showing a change 
in the relative intensities of the Bragg peaks. 



TIP-CHANGING PROCEDURE ON Bi2Sr2CaCu208+5 

Contact between the tip and the Bi-2212 sample can re- 
sults in an unstable tunneling junction and the possible de- 
struction of the approached area, thus decreeing the pok- 
ing procedure used in URu2Si2 and CeCoIns unusable. To 
circumvent this problem we developed a method to change 
the configuration of the tip apex which preserves the stabil- 
ity of the tunneling junction and allows subsequent mea- 
surements on the same area. Scanning over large areas 
(^ 5000 A), it is possible to find small spots (^ 1 A) 
over which the tunneling junction is highly unstable. By 
scanning over them at lower bias set points (^ 20 mV) 
while maintaining a large tunneling current (^ 300 pA), 
the tip sample separation is decreased and the tip is more 
susceptible to changes in its apex configuration. Such a tip 
changing event is characterized by a single spike in the tun- 
neling current. The change in the tip apex configuration is 
confirmed by subsequent topographic measurements. 
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FIG. 2. Average spectra (lower panels) from the dl/dV mea- 
surements used to generate the corresponding On (E) curves in 
Fig. 4(a) of the main text (reproduced in the upper panel). The 
spectra were normalized to the their integral over a range of volt- 
ages (—60 < V < 0). The insets represent the asymmtery param- 
eter of the respective individual topographs acquired simultane- 
ously to the conductance maps. 



EFFECTS OF TIP-CHANGING PROCEDURE ON 
SPECTRA 




Bias [mV] 



Bias [mV] 



FIG. 3. (b) Average spectra from the dl/dV measurements 
used to generate the corresponding On{E) curves in Fig. 4(c) 
of the main text (reproduced in (a)), (c) On{E) measured on 
URu2Si2 below Tho (4.2 K) over the same FOV with the differ- 
ent tips, (d) Average spectra from the dl/dV measurements used 
to generate the corresponding On{E) curves in (c). The insets 
represent the asymmtery parameter of the respective individual 
topographs acquired simultaneously to the conductance maps. 



SETPOINT CONDITION FOR On{E) MEASUREMENTS 

Table I shows the setup current and bias voltage used for 
measuring the conductance maps from which the On{E) 
curves presented in Figs. 3 and 4 of the main text are ob- 
tained. 



The average spectrum for each of the measurements 
in Fig. 4 of the main text is displayed in Figs. 2, 3, and 
4. Notice that the measurements on URu2Si2 displayed 
in Fig. 4(c) of the main text were taken just below Tho 
(15 K), where the effects of the hidden order gap on the 
spectrum are small [4] (see Fig. 3(b)) but where the sensi- 
tivity of On{E) to the hidden order is strong (Fig. 3(a)). 
For comparison we show similar measurements taken at 
4.2 K where the hidden order gap is fully developed. 

Small variations between measurements are expected 
since the different tip-geometries should yield different 
tunneling matrix elements. However, the relative small 
variance over the different measurements indicates that 
tunneling spectra is mostly unchanged by the tip-changing 
procedures used here and that the large variance seen in the 
On (E) curves displayed in Fig. 4 of the main text is likely 
due to different tip geometries. 



TABLE I. 



Material 


Bias (mV) 


Current (pA) 


Location 


CeCoIns 


-200 


1600 


Fig. 3 and Fig. 4(a), cyan 


CeCoIns 


-200 


300 


Fig. 4(a), green 


CeCoIn5 


-200 


300 


Fig. 4(a), red 


Bi-2212 


400 


200 


Fig. 4(b), all 


URu2Si2 


-150 


150 


Fig. 4(c), cyan circles 


URu2Si2 


-150 


150 


Fig. 4(c), cyan inv. triangles 


URu2Si2 


-150 


200 


Fig. 4(c), orange 


URU2Si2 


-200 


200 


Fig. 4(c), red 


URu2Si2 


100 


200 


Fig. 4(c), green 



DRIFT CORRECTION 

All the conductance maps were corrected for drift using 
the analysis method outlined in Ref. [5]. This allowed the 



Bi2Sr2CaCu208^^ 




-100 100 

Bias [mV] 




Bragg peaks in the DFTs of the conductance maps to be 
represented by a single pixel. However, we find that the re- 
sults presented in the current work are not dependent on the 
use of this analysis tool, and the same effects are present 
with the same magnitudes in the raw data. 
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FIG. 4. Average spectra (lower panels) from the dl/dV mea- 
surements used to generate the corresponding On (E) curves in 
Fig. 4(b) of the main text (reproduced in the upper panel). The 
spectra were normalized to the their integral over a range of volt- 
ages (—100 < y < 0). The insets represent the asymmtery pa- 
rameter of the respective individual topographs acquired simulta- 
neously to the conductance maps. 



